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The nondegenerate two-orbital Hubbard model is studied within the dynamic mean-field theory 
to reveal the influence of two important factors, i.e. crystal field splitting and interorbital hopping, 
on orbital selective Mott transition (OSMT) and realistic compound Ca2- a; Sr :E Ru04. A distinctive 
feature of the optical conductivity of the two nondegenerate bands is found in OSMT phase, where 
the metallic character of the wide band is indicated by a nonzero Drude peak, while the insulating 
narrow band has its Drude peak drop to zero in the mean time. We also find that the OSMT regime 
expands profoundly with the increase of interorbital hopping integrals. On the contrary, it is shown 
that large and negative level splitting of the two orbitals diminishes the OSMT regime completely. 
Applying the present findings to compound Ca2- a: Sr a ;Ru04, we demonstrate that in the doping 
region from x = 0.2 to 2.0, the negative level splitting is unfavorable to the OSMT phase. 
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I. INTRODUCTION 

The orbital degree of freedom is one of the key points to 
understand the rich and complicated properties of transi- 
tion metal oxides^, and its roles in Mott-Hubbard metal- 
insulator transition (MIT) have attracted much interest 
in recent years. The orbital selective Mott transition 
(OSMT) in multi-orbital correlated systems, where the 
Mott-Hubbard metal-insulator transition (MIT) in dif- 
ferent orbitals happens successively with the variation 
of doping, interaction, or field, was proposed— to inter- 
pret the two consecutive phase transitions and spin state 
transition in the layered perovskite Ca2- 2: Sr 2 ;Ru04 3 . Be- 
sides, it was also applied to other orbital compounds, 
such as V2O3 and VO2 4 . However, this OSMT scenario 
has not been verified by experiments in the OSMT can- 
didate Ca2- :c Sr. r Ru04 till now. The optical spectra- and 
the angle-resolved photoelectron spectroscopy (ARPES) 6 
experiments were performed to examine the electronic 
states and spectral properties of Ca2- x Sr a; Ru04, but nei- 
ther considerable contraction of Fermi surface volume 5 
nor the OSMT gap for narrow d xz / yz bands^ was found 
near x — 0.5. Whereas the most recent result ob- 
tained by ARPES experiment on Ca2- :E Sr :E Ru04 at x = 
0.2 presented a different scenario, where the low-energy 
quasiparticle excitations in wide d xy orbital disappeared 
completely^. 

Great theoretical efforts have been paid to un- 
earth the essential of the OSMT by the dynam- 
ical mean-field theory (DMFT) and some other 
approache o 8 i 9 i 10 ' 1:L i 12 ' 13 i 14 ' 15 . It has been confirmed that 
the two successive MITs and the OSMT are robust fea- 
tures in the two-orbital and three-orbital Hubbard mod- 
els with asymmetric bandwidths. Meanwhile, the essen- 
tial aspects of the OSMT are not well understood. For ex- 
ample, the effect of Hund's coupling in OSMT is still not 
clear. Koga et al, s and Liebsch et al£ found that the full 
Hund's coupling is crucial for the existence of OSMT in 



two-orbital Hubbard model with asymmetric bandwidth 
using the DMFT with quantum Monte Carlo solver and 
exact diagonalization solver, respectively. However, uti- 
lizing the slave spin mean-field theory and DMFT, Medici 
et a/— - suggested that the OSMT exists over a wider 
range of J when the ratio of the narrow and wide band- 
widths is small enough, even in the full SU(4) spin-orbital 
symmetric case with J = 0. On the other hand, Lieb- 
sch et alJ£ recently introduced a three-orbital Hubbard 
model with positive crystal field splitting for the real 
compound Ca2_ x Sr I Ru04, and found a metallic-band 
insulator transition for the completely filled d xy orbital 
and a Mott transition in the remaining half-filled d xz / yz 
bands, precluding the possibility of the OSMT phase in 
Ca2- x Sr a; Ru04. 

In order to get a clear picture about the unconventional 
MITs in compound Ca2- a: Sr a ;Ru04, all major essentials 
of the real compound should be taken into consideration 
when a realistic model is set up. In layered perovskite 
Ca2- :r Sr :E Ru04, the nearest-neighbor (NN) interorbital 
hopping integrals are comparable with that of the in- 
traorbital integrals in the order of magnitude. However, 
in previous studies^^^ 1 ^ 3 -! 1 ^ 5 -^, the NN interor- 
bital hopping integrals were neglected so as to make the 
self-energy and Green's function be diagonal in orbital 
space. On the other hand, only the hybridization of 
the two orbitals at the same site was considered recently 
in reference [11,17]. It is expected that the presence of 
the interorbital hopping will profoundly affect OSMT. 
Moreover, both the experiments-^ and the first princi- 
ple calculations 1 - demonstrated that with the increase of 
doping in Ca2- K Sr :r Ru04, the structure distortions aris- 
ing from the rotation and tilting of RuO@ octahedra play 
an important role in the two successive MITs. As we will 
analyze later, such distortions lead to the decline of the 
level center of the narrow band with respect to the wide 
one, giving rise to a negative crystal field splitting. 

The aim of this paper is to elucidate the combining 
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effect of two important factors, i.e. interorbital hopping 
and crystal field splitting, on the optical conductivity and 
the OSMT in multiorbital correlated model and realistic 
multiorbital material Ca2_ x Sr x Ru04. We study the opti- 
cal conduction of two-orbital Hubbard model within the 
framework of the extended linearized dynamical mean 
field theory (ELDMFT)^, which has been proved exten- 
sively to be an efficient method, especially for the compli- 
cated multiorbital models. We demonstrate how the in- 
terorbital hopping and the negative crystal field splitting 
affect the nature of OSMT and the optical conductivity of 
the systems. We find that the OSMT regime enlarges in 
the presence of the interorbital hopping, but the negative 
crystal field splitting strongly restrains the occurrence of 
OSMT. In addition, the unconventional MITs in com- 
pound Ca2- x Sr x Ru04 is also discussed. The rest of the 
paper is organized as follows. In Sec. [Til we describe the 
model Hamiltonian and outline the ELDMFT method; 
in Sec. IIII1 we present the major theoretical results of 
the interorbital hopping and the crystal field splitting on 
the OSMT, and the possible application on compound 
Ca2- x Sr :!; Ru04. Sec. [TV]is devoted to the conclusions. 



II. MODEL AND METHODS 

We start from the two-orbital Hubbard model, 
H = Y, E twCl,C iVa + - n)Cf ll7 C iU , 

(ij) IV a ila 
ila ,a 
i,l^l f i,ly£V ,aa' 

(i) 

Here index I (V) = a or b refers to the narrow or wide 
orbital, and t aa (fob) and t a b denote the NN intraorbital 
and interorbital hopping integrals, respectively; U and 
U' represent the on-site intraorbital and the interorbital 
Coulomb repulsions between electrons; and J denotes the 
Hund's coupling. The crystal field splitting arising from 
the structure distortions is defined as = e a — e/,. 

In the presence of the finite t a b, we introduce a canon- 
ical transformation^, 

C la a = ua la + V(3ia 

Ciba = ~va ia + uf3i a . 

Where u and v are the coefficients of the canonical trans- 
formation; a and (3 represent fermion annihilation op- 
erators of the new effective orbitals, with only intraor- 
bital hoppings t a = t aa u 2 + fobV 2 — t a buv and tp — 
taaU 2 +fobV 2 + t a buv , and vanishing interorbital hoppings. 
As a result, the two-orbital model with interorbital hop- 
ping (Eq. (jTJ)) is transformed into a new effective two- 



orbital model without interorbital hopping, 

+ ( £ a - M) J2(at a i<y + P&Pia) 
ia 

~ e d'^(v 2 <x[a a i<y - uv(af a P ia + (3t a a irT ) + u 2 /3+/3 ia ) 

ia 

+H?l f - (2) 

Here Hf-J is the interaction Hamiltonian under the 
canonical transformation. From Eq.(2), we can distin- 
guish the effects of the interorbital hopping from the 
crystal-field splitting more clearly within the new fermion 
representation. That is, the finite interorbital hopping 
tab increases the bandwidth ratio of the two effective or- 
bitals; while the crystal field splitting induces a certain 
extent of the orbital hybridization of the two effective or- 
bitals at the same site. When 6^=0, the Green's function 
is diagonal in the new fermion representations, and we 
can employ the DMFT approach to study the two-orbital 
model easily. On the other hand, when both the crys- 
tal field splitting and the interorbital hopping are taken 
into account, the off-diagonal components of the Green's 
function appear in this fermion representation. In this 
situation, to discuss the effect of crystal field splitting, 
we can use the similar technique adopted by Medici et 
aZjii and Bucnemann et al^ for the two-orbital models 
with the on-site orbital hybridization. In addition, fur- 
ther transformation can be introduced to diagonal the 
Green's function. More details of our method are given 
in Sec. IhTCI 

In the DMFT approach, the self-consistent loop is built 
to fulfil two equations, 

G\ mp (y>) = G\u,) (3) 

and 

S| mp ( W ) = E , ( W ). (4) 

Where G\ mp (w) and T, l imp (uj) represent respectively the 
Green's function and the self-energy of the Anderson 
impurity model mapped from the two-orbital Hubbard 
model. We adopt the ELDMFT approach—., which is a 
fast and efficient method for the multi-orbital systems. 
The result of ELDMFT for one-band Hubbard model 
has been compared with that obtained by the exact 
diagonalizatior*2i and numerical renormalization group 22 
solvers, and a good agreement has been achieved in a 
rather wide U range. 

In the DMFT approach with the exact diagonaliza- 
tion and numerical renormalization group solvers, it has 
been proved that the lower and higher energy behaviors 
of the self-energy play the crucial roles. In order to facili- 
tate the calculation and to consider the lower and higher 
frequency features properly, the ELDMFT approach in- 
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traduces two self-consistent equation; 
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where A fcm = ^ ^ ^+xYk m ( e h <f>j 



The contribu- 
tion of the crystal field Hamiltonian, Hc.f, reads 



Where Vi is the hybridizing parameter in the Anderson 
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bridizmg parameter m tne Anderson r 

1/(1 - dT, l {Q)/duj) is the quaiparticle H CF = drij+(r)V(r)iJj a (r)dr. (10) 



weight, and n (n| mp ) represents the occupation number 
of the lattice (impurity) model. In addition, the free- 
particle DOS is a semicircular type. These two equations 
correspond to those of the conventional DMFT approach 
in Eq. ([4]), respectively, and there is no other extra re- 
striction employed in ELDMFT. The ELDMFT method 
is able to access the entire parameter region since the 
Luttinger theorem (N) — nL (L is the number of lattice 
size) is exactly fulfilled in the whole range of filling. 

In the new fermion representation obtained by the 
canonical transformation 2 ^, the optical conductivity of 
each effective orbital in the paramagnetic system with- 
out the interorbital hopping read s 21 ' 23 
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= 7T 



duj'deD l {e)p l e {u)')p l e {uj' +u) 



(7) 



where D l (e) and p l (cu) (p l (uj) = —^ImG l (uj + ij)) de- 
note the transition matrix elements and the density of 
states (DOS) in the new orbital channel a or f3, respec- 
tively; and f(u>) is the Fermi-Dirac distribution function. 
A small spectra broadening (7=0.1) is introduced for the 
calculations of DOS. When there is no hopping between 
the new quasiparticle channels, the total optical conduc- 
tivity <t(ui) is the summation of the contributions of two 
quasiparticle orbitals, a a (u>) and 073 (w). In the presence 
of interorbital hopping, the full optical conductivity for 
the Green's functions and the transition matrix D u (e) 
with off-diagonal terms can be obtained by, 

a(w) = tt f dw'deY^D 11 ' {e)p{{u')p{{J + u) 
^ l,l> 

x /("')-/("' + "). (8) 
u> 

In this paper, we mainly focus on the optical conductiv- 
ity of each orbital, 07, so as to elucidate the metallic or 
insulating nature of each orbital. 

In many transition-metal oxides, the transition-metal 
ions and oxygen ions usually form TM0 6 octahedra. 
The TM06 clusters distorted from the perfect octahe- 
dra contribute the crystal field splitting. Starting from 
the single-particle potential V(r), one can estimate the 
splittings in the point-charge approximation of the crys- 
talline field theory 2 ^. Expanding the single-particle po- 
tential V(r) = J^j — ^1 in terms of the spherical 



In the second quantization representation, Hcf can be 
rewritten as 

j ll'a 

where 



V w = / w l Jr)V(r)wUr)dr, 



(12) 



with ipcr(r) = J2iJ2j w j(r) c ]i<? an d Wj being the Wan- 
nier basis consists of the wavefunctions of the transition- 
metal d orbits and of the nearest-neighbor atomic orbits. 
In the high symmetric situation, the matrix Vw is diag- 
onal, giving rise to e/ = Vw . Later, we will evaluate the 
crystal field splitting based on these equations and the 
experimental data of the crystal structures of the family 
compounds Ca2- 2; Sr 2 ;Ru04. 



III. RESULTS 

We study the effects of interorbital hopping and crystal 
field splitting separately in the asymmetric model, and 
discuss the application to real compound Ca2_ a; Sr 2 ;Ru04. 
The calculations are performed for the Bethe lattice with 
semicircular DOS, and our results are measured in units 
of the hopping integral t^b ■ In our calculation, the band- 
widthes of the free-particle spectra are fixed as W a — 2 
and W b = 4. 



A. OSMT and optical conductivity 

To present the distinction of the optical conductivity 
in different phases, we first investigate the simplest two- 
orbital model without interorbital hopping (t a b = 0) and 
crystal field splitting (e^ = 0). The quasiparticle weight 
Z is a powerful tool to characterize the MIT and the 
OSMT. We present the evolution of quasiparticle weight 
with the interaction U in Fig.[lja) for a two-orbital model 
with different orbital bandwidths W1/W2 — 0.5. We find 
that when the Coulomb interaction U increases to U Cl — 
1.45, the first MIT, i.e. the OSMT, happens. The second 
transition occurring at U C2 — 2.90 is the conventional 
MIT. 



4 



N .5 





U' 



UP 



FIG. 1: (Color online) (a) Quasiparticle weight as a function 
of interaction U; (b)-(d) evolutions of optical conductivity in 
different orbital channel with U from metallic (b), OSMT (c), 
to insulating (d) phases in the two-orbital half-filling Hubbard 
model with J = 0.5, taa = 0.5, tbb = 1.0, t a b = 0, and e<j=0. 



FIG. 2: (Color online) Phase diagram of the two-orbital 
Hubbard model (a) without interorbital hopping, and (b) 
with interorbital hopping. Theoretical parameters: J — 0.5, 
t aa = 0.5, t bb = 1.0, and e d =0. 



The quasiparticle contribution to the optical conduc- 
tivity arises from the electronic transitions between the 
lower and the upper Hubbard subbands. As shown in 
Fig. QJb) , the Drude peaks of both orbitals are consider- 
ably large as U < U C1 , which demonstrates that the sys- 
tem is in the metallic phase. On the contrary, the Drude 
peaks of both ortitals vanish completely in the insulating 
region with U > U C2 , as shown in Fig. [Tfd) . However, 
in the OMST regime with U Cl < U < U C2 , the Drude 
peak of the narrow band vanishes and only the wide band 
contributes to the conduction channel in the optical con- 
ductivity, as shown in Fig. [He). Similar to the orbital 
ordering in manganites and vanadates, the OSMT phase 
is also orbital polarized in the wide channel^. Therefore, 
the OSMT phase can be detected by using the polarized 
synchrotron radiation light to measure the distinct opti- 
cal conduction behaviors of different orbitals. 

Apart from the low-energy Drude-like feature coming 
from the coherent quasiparticle DOS existing near Ep, 
there is another spectral feature, as shown in Fig. 1(b)- 
(d), centered at the intermediate energy, which repre- 
sents the charge transfer from the lower and upper Hub- 
bard subbands. Since the general feature of DOS, i.e. 
the lower and upper Hubbard subbands and a quasi- 
particle resonance at Fermi surface, can be reproduced 
by ELDMFT, our numerical results of optical conductiv- 
ity are qualitatively accurate. To quantitatively compare 
with experiments on Ca2-xSr K Ru04, more accurate DOS 
should be adopted rather than the semicircular DOS, and 
all three t2 g orbitals should be taken into account. This 
needs to be studied in future work. 



B. Effect of interorbital hopping 

Since the interorbital hopping is significant in real ma- 
terial, it is crucial to dig out its effect on the phase di- 
agram and the optical conductivity. To our knowledge, 



no work has been done to take this important factor into 
consideration. Till now, only the hybridization of the two 
orbitals at the same site had been studied by DMFT 17 , 
where the OSMT tends to be destroyed. However, the 
physical origins of the on-site intra-atomic hybridization 
and NN interorbital hopping are quite different. By an- 
alyzing the quasiparticle weight at Fermi surface, we ob- 
tain the phase diagram of the two-orbital Hubbard model 
with interorbital hopping integrals t a \, = 0.5, as shown in 
Fig. HJb). For the convenience of comparison, the phase 
diagram for t a b = 0.0 is also shown in Fig. [2a). In both 
cases, a stable OSMT regime (the red shadow) lies be- 
tween the Mott insulating (denoted by " I" ) and the con- 
ventional metallic (denoted by "M") phases. When the 
interorbital hopping is finite, for example t a t = 0.5, the 
area of the OSMT regime greatly expands, in comparison 
with the case t a h = 0. 

It should be noted that the OSMT we observed here 
is for the effective orbitals in the new fermion represen- 
tation, which are the hybridized orbitals of the two re- 
alistic orbitals^. Within this representation, one could 
find that the bandwidth ratio of the current two effec- 
tive orbitals increases with the increase of the interor- 
bital hopping integrals t a \,\ the wide quasiparticle band 
becomes wider, while the narrow one becomes narrower. 
Thus the OSMT occurs more easily, and its area expands 
a lot in the phase diagram for large t a b- In addition, we 
also find that for zero and nonzero interorbital hopping, 
the OSMT regime is close to the U = V line, but never 
crosses it. Our result of t a {,=0 agrees with Medici et al.'s 
phase diagram obtained by slave boson approach; how- 
ever, at J — 0.5, our result is slightly different from the 
result obtained by Koga et al. 8 for J = in which the 
OSMT regime crosses the U = U' line. 

The important role of the interorbital hopping t a b in 
the OSMT also manifests in the optical conductivity. 
Starting from an insulating case with U — 2.25, U' = 
and t a b = 0, we obtain the typical feature of insulating 
optical conductivity for both the narrow band and the 
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FIG. 3: (Color online) The optical conductivities of the nar- 
row orbital (a) and of the wide orbital (b) as functions of 
energy for different interorbital hopping integrals: t a b=Q, 0.5, 
and 1.0 as U = 2.25; (c) Quasiparticle weight as a function of 
J with tat = t a - Theoretical parameters: U' = 0, J = U/2, 
t aa = 0.5, t b b = 1.0, and e d =0. 



wide one. Without the interorbital hopping, the system 
is a conventional Mott insulator, and the Drude peaks 
of both effective orbitals vanish. As shown by the black 
dashed lines in Fig.[3{a) and[3fb), there is no Drude peak 
for both <j a (oj) and 073(0;), and the peaks in the interme- 
diate energy region clearly show the hopping between 
the subbands separated by the intraorbital Coulomb in- 
teractions. It is obvious that the interorbital hopping 
t a b does not significantly change the intermediate-energy 
band structures in the DOS of the narrow band. How- 
ever, in Fig.[3f b), we find that the Drude peak of the opti- 
cal conductivity of the wide band 073(0;) crucially changes 
with the increase of the interorbital hopping t a \,. As t a b 
increases to 0.5 and unity, the Drude peak in 073(0;) con- 
siderably increases, implying that the system enters into 
the OSMT regime. We also find that the quasiparticle 
weight drops abruptly at the two successive MIT points, 
as shown in Fig.[3]Jc), when the finite interorbital hopping 
is taken into account. We observe a small plateau near 
the metal-OSMT transition point, which primarily shows 
that both of the two successive MITs are first-order. The 
order of the OSMT deserves further study. 



C. Effect of crystal field splitting 

The level splitting of the two orbitals e d arising from 
the crystal field is another important factor that influ- 
ences the essential of the MIT and the OSMT in the 
two-orbital systems. To study the evolution of the optical 
conductivity of the compound Ca2_ a; Sr a ;Ru04 with the 
crystal field splitting e d , we start from an OSMT phase 
with the parameters U — 2.0, U' = 1.0 and J = 0.5, 
which lies in the red shadow regime in the phase diagram 
shown in Fig.^b). 

As shown in Eq. @, there are hybridization terms 



ufaPia and fifoOUc, for the two effective orbitals on the 
same lattice when both the interorbital hopping and crys- 
tal field splitting are nonzero. In this condition, the 
Green's function is defined as, 



Go (u,k)~ l = 



oj + (J,— e a (k) 

C d UV UJ 



e d uv 



+ fi-e (k) 



, (13) 



where e a (k) and e,a(fc) are the bare dispersion relations 
for the a and (3 effective orbitals, respectively. Obviously, 
we can directly define the local Green function as 



G a{p \uj) = 



d.r 



oj + n - x - E<*W (w) - A°(ffl (w, x) ' 

(14) 



with 



A Q(/3) (u,x) 



( 



e d uv 



+ n - x - (uj) ' 



(15) 



Here £ Q (^(o;) and E q/3 (oj) represent the diagonal and 
off-diagonal terms of the self-energy matrix, respectively. 

Based on the special feature of the above off-diagonal 
Green's function in Eq. (|13[) . we can introduce an unitary 
matrix U(u>, k) to diagonalize Eq. (|13p as an effective 
Green's function. 



Go^r 1 = c/Go^fcr 1 ^ 1 

oj + /i-e+(fc) 

uj + [i — (-(k) 



(16) 



with 



e ± (k) = l(e a (k)+e p (k)±y/^k) - e {k)f + Aeju^). 

(17) 

For a half-filling asymmetric two-orbital system, when 
the relationship \e a (k) — ep{k)\ « |e^| is satisfied in the 
energy region close to the Fermi surface, we can approx- 
imate the above dispersion relations in the first order as, 



W = \[t a {k)+e p {k) 



±(2\e d uv\ 



>(k)~ep(k)f 
<i\e d uv\ 



)]■ (18) 



It is obvious that the finite crystal field splitting mixes 
the two effective a- and /3-orbital in the low energy 
regime, and as a result, the OSMT is suppressed. 

We employ the ELDMFT approach 19 to numerically 
solve the corresponding impurity problem, and find that 
the OSMT phase is strongly influenced by the transfer 
of electrons between the two effective orbitals. We start 
from an OSMT phase with e d = 0.0, in which the elec- 
trons in the narrow orbital are localized and itinerant in 
the wide one, respectively. As shown in Fig. |Uc), l ne 
OSMT character clearly manifests in the Drude peak of 
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FIG. 4: (Color online) Effect of crystal field splitting on op- 
tical conductivity of each orbital and OSMT phase in two- 
orbital Hubbard model with U = 2.0, U' = 1.0, J = 0.5, 
t aa = 0.5, ttb = 1.0, and tab = 0.5. The black solid line de- 
notes the effective narrow orbital, and the red dash line the 
effective wide obital. 



each orbital: the Drude peak vanishes completely for in- 
sulating narrow orbital, while the wide orbital is metallic 
with a finite Drude peak. 

When the crystal field splitting becomes negative as 
e d = -0.5 (see Fig. Ha)) and -0.25 (see Fig. gfb)), the 
narrow orbital shifts to the lower energy region, and a 
fraction of electrons in the wide orbital transfers to the 
narrow orbital. The electron filling in the narrow or- 
bital is more than a half, thus the narrow orbital also 
becomes metallic. In this situation, the Drude peak of 
the narrow orbital appears again. On the other hand, as 
the crystal field splitting is positive, the narrow orbital 
keeps its insulating character and the wide one is metal- 
lic. So the system remains in the OSMT phase, as shown 
in Fig.Htd). 

Furthermore, we find that in both the metallic and 
insulating phases near the border of the OSMT regime, 
a small negative crystal field splitting can drive neither 
the metallic system nor the insulating system into the 
OSMT phase. Therefore, the negative crystal field split- 
ting is not in favor of the stable OSMT phase, even in the 
presence of the finite interorbital hopping integrals. This 
result is contrary to that of Dai et alM- , who supported 
that the OSMT exists in the three-orbital Hubbard model 
with crystal field splitting. 



D. Application to Ca2_ a; Sr a ;Ru04 

For the compound Ca2- a; Sr 2 ;Ru04, the crystal field 
splitting arises from the structure distortions, and it can 
be approximately modeled as a two-orbital system since 
the double degeneracy of the narrow d xz / yz bands is kept 
over the whole doping range. Within the framework of 
the single-ion crystal field theory 2 ^, we investigate the 
crystal field splitting between the d xz / yz orbitals and the 
d xy orbital with the doping concentration x varying from 
0.0 to 2.0. According to the crystal structures for x=2.0, 
0.5, 0.2, 0.1 and 02£, we find that in the most doping re- 
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FIG. 5: (Color online) (a) Crystal field splitting as a 
function of x in compound Ca2- :E Sr I Ru04; (b)Quasiparticle 
weight as a function of td- Other parameters are the same to 
Figure 4. 



gion, the narrow d xz / yz bands are lower than the broad 
d xy band, showing that the negative crystal field splitting 
is dominant. Noting that though the single- ion crystal 
field theory is relatively simple, it is qualitatively correct 
for the sign and the relative magnitude of the crystal 
field splittings. As shown in Fig. 02a), the averaged split- 
tings between the d xz / yz and the d xy bands are ed/td = 
-0.013, -0.11, -0.012, 0.0025 and 0.005 for x=2.0, 0.5, 0.2, 
0.1 and 0, respectively. Here ed is a positive constant, 
representing the 4ti-orbital radial integral of the single- 
particle potential (e d[) = r 2 dr(4>(r)\V (r)\4>(r)}). The 
level separation of the narrow d xz / yz bands from the wide 
d xy band in Ca2- a ;Sr 2 ;Ru04 varies gradually from posi- 
tive to negative with the increase of Sr concentration, and 
have the largest value as x — 0.5. This largest level sep- 
aration arises from the rotation- type distortion 2 ^, which 
becomes the most significant when x = 0.5. As x further 
decreases, the tilting distortion of the RuO@ octahedra 26 
becomes more important, which usually mixes different 
orbitals. 

Since the normal state of Sr2Ru04 is metallic, our pre- 
ceding results shows that negative crystal field splitting 
is unfavor of stable OSMT phase. Therefore, our re- 
sults exclude the possibility of the OSMT occurring in 
Ca 2 - 2 ;Sr 2 ;Ru04 when x increases from 0.2 to 2.0. Very 
recently, Liebsch et al^ suggested that in compound 



Ca2- a; Sr a ;Ru04, the splitting of the d x 



</yz lo the d xy 
bands is positive. In their study, the band structures in 
compound Ca2- a; Sr a ;Ru04 is rigid, and thus the variation 
of electron occupations in different orbitals arises from 
the change of the positive crystal field splitting. Based 
on this hypothesis, Liebsch et al~ and Werner et al£Z- 
discussed the role of the positive crystal field splitting on 
the MIT in two-orbital Hubbard model. However, Lee et 
al£ have demonstrated that with the doping increasing 
from x=0 to 2, the bandwidths of the d xy and d xz / yz or- 
bitals considerably change, implying that the rigid band 
approximation is not proper. 



7 



Based on the above findings, the two successive MIT in 
Ca 2 _2;Si' :I ;Ru04 is obviously independent of the OSMT. A 
question arises: what happens in the first unconventional 
MIT at x=0.2? We notice that at the doping concentra- 
tion of x=0 and x=0.5, the doped ruthenates are found 
to be orbital ordering. Therefore, the two successive MIT 
should be accompanied by the transitions of the ordered 
orbital states, which have not been well treated in the lit- 
erature. To account for the properties of these two MIT, 
the orbital ordering should be taken into account prop- 
erly within the DMFT approach, which has gone beyond 
the scope of this paper. 

IV. CONCLUSIONS 

We have studied the effects of two important factors, 
i.e. the crystal field splitting and interorbital hopping, on 
OSMT in nondegenerate two-orbital models and in the 
compounds Ca2_ a; Sr :I ;Ru04. Our results show that the 
negative crystal field splitting can hybridize the two non- 



degenerate orbitals at the same lattice site, and strongly 
restrain the appearance of the OSMT phase. However, 
the finite interorbital hopping considerably expands the 
OSMT region in the phase diagram. In the doping region 
form x = 0.2 to 2.0 in C^-zSr^RuO-i, the crystal field 
splitting introduced by the distortion of the RuOg octa- 
hedra is found to be negative. Therefore, the occurrence 
of OSMT in compound Ca2- a; Sr a ;Ru04 unlikely happens 
in the corresponding doping region. 
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